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Abstract 

In analogy to gas-dynamical detonation waves, which consist of a 
shock with an attached exothermic reaction zone, we consider herein 
nonlinear traveling wave solutions, termed "jamitons," to the hyper- 
bolic ("inviscid") continuum traffic equations. Generic existence cri- 
teria are examined in the context of the Lax entropy conditions. Our 
analysis naturally precludes traveling wave solutions for which the 
shocks travel downstream more rapidly than individual vehicles. Con- 
sistent with recent experimental observations from a periodic roadway 
(Sugiyama et al. New Journal of Physics, 10, 2008), our numerical cal- 
culations show that, under appropriate road conditions, jamitons are 
attracting solutions, with the time evolution of the system converging 
towards a jamiton-dominated configuration. Jamitons are character- 
ized by a sharp increase in density over a relatively compact section 
of the roadway. Applications of our analysis to traffic modeling and 
control are examined by way of a detailed example. 

PACS: 89.40.Bb Land transportation; 47.10. ab Conservation laws; 
47.40.Rs Detonation 

1 Introduction and problem formulation 

The economic costs in terms of lost productivity, atmospheric pollution and 
vehicular collisions associated with traffic jams are substantial both in de- 
veloped and developing nations. As such, the discipline of traffic science has 



expanded significantly in recent decades, particularly from the point of view 
of theoretical modeling [1]. Borrowing terminology applied in Payne [2] and 
elsewhere, three generic categories describe the approaches considered in 
most previous analyses. "Microscopic" models, such as "follow the leader" 
studies [3] or "optimal velocity" studies [4] consider the individual (i.e. La- 
grangian) response of a driver to his or her neighbors, in particular, the 
vehicle immediately ahead. "Mesoscopic" or "gas-kinetic macroscopic" anal- 
yses, such as the examinations of Phillips [5] and Helbing [6] take a statistical 
mechanics approach in which vehicle interactions are modeled using ideas 
familiar from kinetic theory. Finally "macroscopic" studies [2, 7-12] model 
traffic flow using conservation laws and a suitable adaptation of the methods 
of continuum mechanics [13, 14], which yields governing equations similar to 
those from fluid mechanics. It is this latter category of analysis that is of 
interest here. 

Treating the traffic flow as a continuum, we begin by considering a one 
dimensional Payne- Whitham model with periodic boundary conditions, i.e. 
vehicles on a circular track of length < A < oo [15]. The governing 
equations for mass and momentum are then (e.g. Kerner &; Konhauser [10]) 



where the subscripts indicate differentiation, r is a relaxation time-scale, u is 
the traffic speed, and p is the traffic density — with units of vehicles/length. 
The traffic pressure, p, which incorporates the effects of the "preventive" 
driving needed to compensate for the time delay r, is typically assumed 
to be an increasing function of the density only, i.e. p = p{p) [11]. Here, 
in order to have a well-behaved theoretical formulation in the presence of 
shock waves [16], wc shall assume that p is a convex function of the specific 
volume V = 1/p (road length per vehicle). This implies that dp/dv < 
and d^p/dv^ > 0, which holds for the functions typically assigned to p in 
macroscopic models. Finally, u = u{p) gives, for a particular traffic density, 
the desired or equilibrium speed to which the drivers try to adjust. The 
precise functional form of u is, to a certain degree, rather arbitrary and 
indeed several variants have been proposed [1]. Generally, tt is a decreasing 
one-to-one function of the density, with < u{0) = uq < oo and u{pm) = 0, 
where: 

• pm denotes the maximum density, at which the vehicles are nearly 
"bumper-to-bumper" — thus i = p'^ is the "effective" (uniform) ve- 
hicle length. 
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• uq is the drivers' desired speed of travel on an otherwise empty road. 

In this paper, we restrict ourselves to the representative form u = uq{1 — 
p/pm)^, where n is "close" to 1. We defer the detailed treatment of the 
exact conditions on u that guarantee the existence of self-sustained nonlinear 
traveling waves in traffic (termed "jamitons" herein) to a later publication. 

Ubiquitous attributes of the solutions to continuum traffic models are 
stable shock-like features (see e.g. Kerner & Konhauser [9] and Aw & Ras- 
cle [11]). In analyzing such structures, a dissipative term proportional to 
Uxx, analogous to the viscous term in the Navier-Stokes equations, is often 



added to the right-hand side of the momentum equation (1.2) in order to 
"smear out" discontinuities [1]. However, the physical rationale for this term 
is ambiguous and the proper functional form is therefore subject to debate. 
Solutions, such as those obtained by Kerner & Konhauser [9], whose dy- 
namics are non-trivially influenced by viscous dissipation must therefore be 
interpreted with care. Herein, an alternative line of inquiry is proposed: we 
seek self-sustained traveling wave solutions to the "inviscid" equations ( |1.1| 



- 1.2 ) on a periodic domain, where shocks are modeled by discontinuities, as 
in the standard theory of shocks for hyperbolic conservations laws [16]. As 
we demonstrate below, not only do such nonlinear traveling waves exist, but 
they have a structure similar to that of the self-sustained detonation waves 
in the Zel'dovich-von Neumann-Doering (ZND) theory [17]. According to 
the ZND description, detonation waves are modeled as shock waves with an 
attached exothermic reaction zone. In a self-sustained detonation wave, the 
flow downstream of the shock is subsonic relative to the shock, but accel- 
erating to become sonic at some distance away from the shock. Hence, the 
flow behind a self-sustained detonation wave can be "transonic," i.e. it may 
undergo a transition from subsonic to supersonic. The existence of the sonic 
point, the location where the flow speed relative to the shock equals the local 
sound speed, is the key feature in the ZND theory that allows one to solve 
for the speed and structure of the detonation wave. Its existence also means 
that the shock wave cannot be influenced by smooth disturbances from the 
flow further downstream so that the shock wave becomes self-sustained and 
independent of external driving mechanisms. Hence the sonic point is an 
"acoustic" information or event horizon [18]. 

For the nonlinear traffic waves to be discussed herein, this means that 
their formation, due to small initial perturbations, is analogous to the ig- 
nition and detonation that can occur in a meta-stable explosive medium. 
Although this analogy has not, to our knowledge, been reported in the 
traffic literature, the physical and mathematical similarities between deto- 
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nation waves and hydraulic jumps, described by equations similar to (1.1) 



and ( 1.2 ) , were recently pointed out by Kasimov [19] . (The analogy between 
hydraulic jumps and inert gas-dynamic shocks was recognized much earlier 
~ see e.g. Gilmore et al. [20] and Stoker [21].) As with Kasimov's analysis, 
our aim is to herein exploit such commonalities to gain additional under- 
standing into the dynamics of traffic flows, in particular, the traffic jams 
that appear in the absence of bottlenecks and for no apparent reason. From 
this vantage point, novel insights are discerned over and above those that 
can be realized from the solution of a Riemann problem [11] or from the lin- 
ear stability analysis of uniform base states (see e.g. Appendix [C]) . Indeed, 
when such linear instabilities are present initially, our extensive numerical 
experiments suggest that the resulting "phantom jams" (see Helbing [1] and 
the many references therein) will ultimately saturate as jamitons. This ob- 
servation provides a critical link between the initial and final states, the 



latter of which can, under select conditions (see e.g. § 5.3), be described 
analytically. Moreover, as we plan to illustrate in forthcoming publications, 
such self-sustained traffic shocks are also expected on non-periodic roads. 
Thus the model results presented below can be readily generalized beyond 
the mathematically convenient case of a closed circuit. 

The rest of the paper is organized as follows: in § [21 we outline the 



basic requirements for ( 1.1 ) and ( 1.2 ) to exhibit traveling wave solutions. To 
demonstrate the generality of this analysis, we consider in §[3] modified forms 
for the momentum equation (developed by Aw & Rascle [11] and Helbing 
[12]). From this different starting point, the salient details of § [2] shall be 
reproduced. The analysis is further generalized in § |4] which considers a 
phase plane investigation of the governing equations from § |2] and § |3] A 
particular example is studied, both theoretically and numerically, in § |5] in 
which u and p are assigned particular functional forms. The impact of our 
findings on safe roadway design is briefiy discussed. Conclusions are drawn 
in §[6| 

2 Traveling wave solutions — jamitons 

To determine periodic traveling wave solutions to the traffic fiow equations 



(1.1) and ( |1.2[ ), we begin by making the solution ansatz, p = p{r]) and 
u = u {rj) , where the self-similar variable rj is defined by 

X — st , , 

r?= . (2.1) 
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Here s is the speed, either positive or negative, of the travehng wave. Equa- 



tion (1.1 ) then reduces to 



m 



(2.2) 



where the constant m denotes the mass flux of vehicles in the wave frame 



of reference. Substituting (2.2) into (^1.2|, we obtain 



du {u — s) {u — u) 
drj 



{u 



^2 2 
s) — & 



(2.3) 



where we interpret u = u{p) as a function of u via (^. Herec= (pp)^/^ > 



is the "sound speed," i.e. the speed at which infinitesimal perturbations move 
relative to the traffic flow. 



Equation (2.3) is a first order ordinary differential equation and there- 



fore, barring pathological and unphysical choices for c and u, does not admit 
any smooth periodic solutions. Hence, the periodic traveling wave(s) — if 
they exist — must consist of monotone solutions to (2.3) that are connected 



by shocks. The simplest situation, as reproduced in the physical experiment 
of Sugiyama et al. [15], is one in which there is exactly one shock (with speed 
s) per period. The case of multiple shocks per period is more complex. We 
briefly address this situation in 



5.3 



Before going into the details of the solution, we recall that shocks must 
satisfy two sets of conditions to be admissible. First, they must satisfy the 
Rankine-Hugoniot conditions [11, 13, 14], which follow from the conservation 
of mass and momentum, and ensure that shocks do not become sources 
or sinks of mass and/or momentum. For the equations in ( |1.1| - [T^ , the 
Rankine-Hugoniot conditions take the form. 



s[p] 
s[pu\ 



[pu] 

[p + pu^\ 



(i.e. conservation of mass), (2-4) 
(i.e. conservation of momentum), (2.5) 



where s is the shock speed and the brackets [ ] indicate the jump in the 
enclosed variable across the shock discontinuity. These equations relate the 
upstream and downstream conditions at the shock. In particular, let the 
superscripts -|- and — denote the states immediately downstream (right) 



and upstream (left) of the shock, respectively. Then (2.4) is equivalent to 
p'^ {u^ — s) = m = p^ {u^ — s), (2-6) 



where the constant m is the mass flux across the shock. Of course, for a 



shock embedded within a jamiton, this m is the same as the one in (2.2) 
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Figure 1: Characteristics on both sides of a left (a) and right (b) shock in 
the frame of the shock. The flow direction is from left-to-right in (a) and 
right-to-left in (b). The equations of the characteristics are C± : drj/dt = 
(li — s lb c)/t, where ry is defined by (2.1 ), and s is the shock speed. 



Second, shocks must satisfy the Lax "entropy" conditions [14, 16], which 
enforce dynamical stability. In the case of shocks in gas-dynamics, these 
conditions are equivalent to the statement that the entropy of a fluid parcel 
increases as it goes through the shock transition — hence the name. How- 
ever, the existence of a physical entropy is not necessary for their formu- 
lation: stability considerations alone suffice. Furthermore, these conditions 
also guarantee that the shock evolution is causal. 

For the particular system of equations in ( |1.1| - 11.2[ ) , the Lax entropy 



conditions — given below in (2.7 - 2.8) — state that one family of charac- 
teristic^ must converge into the shock path in space-time, while the other 
family must pass through it. Thus exactly two families of shocks are possi- 
ble, as illustrated in figure [T| The left (respectively right) shocks have the 
left (respectively right) characteristics converging upon them. 

In the context of the periodic traveling waves with a single shock per 
period, the above discussion implies that, in principle, two cases are possible: 
jamitons containing a left shock wherein the mass flux, m, is positive ~ see 
item [T] below, and jamitons containing a right shock wherein the mass flux, 
m, is negative ~ see item [2] below. However, as we argue after item [2] 
only jamitons with m > are mathematically consistent; self-sustained 
traveling waves carrying within them a right shock are not permitted. This 
is consistent with the experiment of Sugiyama et al. [15], as well as with 
one's everyday driving experience: it is situations where individual vehicles 
overtake shocks (hence m > 0) that are observed in reality, rather than the 
converse. Thus whereas in second order traffic models information in the 
form of shock waves can travel downstream faster than individual vehicles 
[22], the results in this paper show that this cannot happen in the form of 



^The curves in space-time along which infinitesimal perturbations propagate. Namely: 
dx/dt = u + c ("right" characteristics), and dx/dt = u — c ("left" characteristics). 
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a self-sustained traveling wave. Further to the analysis of Aw & Rascle [11] 
and Helbing [12] , this observation lends additional support to the conclusion 
that second order models are not ipso facto flawed — see also Helbing [1], 
§ III.D.7 and the references therein. 

Let us now consider in some detail the two scenarios that can, in princi- 
ple, arise for traveling waves with a single shock per period. 

1. The shock is the left shock, as in figure [T]a, so that 

(u-c)~ > s > (u-c)+. (2.7) 

In this case the mass flux must be positive, since m = [u^ — s) > 
> 0. Moreover, p'^ / p^ = {u — / {u — s)^ > c~ /c^ , so that 
p+c"*" > p~c~ . The traffic pressure p is a convex function of v = p~^ 
(see §[!]), hence pc= (— dp/dv)^/^ is an increasing function of the traf- 
fic density. Thus p^ > p^ implies that p^ > p~ . In other words, 
the shock is compressive: the traffic density increases as the vehicles 
pass through the shock, traveling from left to right. Conversely, since 
u = s + m/ p, it follows that u~ > u'^ and, consequently, vehicles de- 
celerate as they overtake the shock. It should then be clear that these 
shocks have all the familiar properties of traffic jams. 

We conclude that for a traveling wave with a left shock, the continuous 



solution of (2.1 - 2.3) must have a decreasing density, dp/drj < 0, and 
an increasing velocity, du/dr] > 0. This follows because the solution 
must be a monotone function of ry, and must connect the post-shock 
state (yO^, n^) in one shock, to the pre-shock state (p~, u~) in the 
subsequent shock across a period in r/ — say from ry = to r/ = A. 

2. The shock is the right shock, as in figure [T]b, so that 

{u + c)~ > s > (u + c)^ . (2.8) 

The mass flux now is negative, since m = p^ {u — s)~^ < —p~^ c+ < 
0. As in item [ij it is straightforward to show that u" > and 
p~ > p'^. For a traveling wave with a right shock, the continuous 



and monotone solution of (|2.1|- 2.3) must have both the density and 



velocity increasing with rj (i.e. dp/drj > and dn/d?7 > 0), in order 
to connect (p+, m+) to (p~, u~) across a period in ry. Notice that the 
shock is again compressive: the traffic density (in the vehicles' frame 
of reference) increases as vehicles pass through the shock. However, in 
this latter case, the shock overtakes the vehicles from behind, which 
accelerate as they pass through the shock transition. This is a clearly 
counter-intuitive situation, not observed in real traffic [22]. 
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Fortunately, as we demonstrate next, traveling wave solutions with m < 
are mathematically inconsistent, which obviates the need to consider them 



any further. First, (2.2) is employed to rewrite (2.3) in the form 

du 

--- = mpG{p, s, m), (2.9) 
ari 

where 

^ u — u , , 

G=—2 2.10 

— = p^ {(m — s)^ — c^}, and u = s + m/ p. Because du/drj > 0, a 
smooth solution connecting (/j"*", u'^) to (p^, u^) requires G < 0. However, 
m? — p^ (? = p{u — s + c) (m — pc) and m — pc < 0, and it follows from (2.8 ) 



that (m^ — c^)+ > > (m^ — p^ c^)~. Thus G < requires {u — u)~^ < < 
{u — u)~ , which is impossible since ■u(p) — u = u{p) — s — m/p is a strictly 
decreasing function of p: u{p) decreases with increasing p by assumption 
and m < 0. 

The difficulties documented in the previous paragraph are avoided for 



traveling waves with m > 0. In this case, we demand a solution of (2.9) 
with du/drj > and dp/dr/ < 0, connecting (p^, it"*") to (p~, u~). This 
in turn requires G > for p~ < p < p^. The assumptions on p = p{p) 
imply that — p^ is a strictly decreasing function of p (item [l| . Since 
rm? — p^ (? = p(u — s — c) (m + pc) with m + pc > 0, it follows from (2.7) that 



(m^ — p^ c^) > > (m^ — p^ c^)"*". We conclude therefore that — (? (? 
has a unique, multiplicity one, zero at some ps, with p^ < ps < p^ ■ In 



order for (2.9) to have a smooth solution with the desired properties, the 



numerator u — u = u{p) — s — mj p must be such that 

(a) u — u has a simple zero at p = ps] 

(b) u-u>0 for p^ < p < Ps] } (2-11) 

(c) u — u < for Ps < P ^ P^ ■ 



This then guarantees not only that the zero in the denominator of (2.9 



2.10) is cancelled by a zero of the same order in the numerator, but that 
the resulting regularized ordinary differential equation yields du/dr] > 
everywhere, as needed. Indeed, ps is the traffic density at the sonic point, 
with corresponding flow speed Us = s + m/ps. 

Figure [2] illustrates the situation, with plots of u and u as functions of p 
and u for representative initial conditions p"*" and u^. The shock speed, s, is 
here restricted by the inequalities n"*" > s > {u — c)^. Clearly, the conditions 
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(a) m>0 



(b) m>0 



^ 0.5 





"Aj/So 0.5 1 
u/uo 

Figure 2: Panel (a) shows typical profiles for u/uq (solid curve) and u/uq 
(dashed curve) as functions of p/pM, where u = s + m/p. Equivalently, 
these profiles may be plotted against u/uq by employing the transformation 
(2.2), as shown in panel (b) — the dashed curve is then just a line with 
unit slope. The two curves will either not intersect at all, be tangent at a 
single point, or have two transversal intersections. The case of interest to 
us is the one with two transversal intersections, as depicted here: U2 is the 
intersection with a smaller velocity and a larger density p2 = m/{u2 — s) — 
solid circles, while the other intersection (open circles) defines ui and pi. It 



should be clear that, in order to satisfy the conditions in (2.11), the sonic 
point must coincide with (^2, P2)- Finally, a physically meaningful solution 
requires u > everywhere. Hence, n2 > is needed. As shown by panel (b), 
this condition is equivalent to the statement: when tt = 0, that \s p = pM, 
the corresponding u = um = s + m/pM is positive. 



in (2.11) require that the sonic point values of the density, ps, and velocity. 



Us, coincide with p2 and U2- The sonic condition therefore reads 

U2 = s + c{p2), 



(2.12) 



from which the jamiton speed, s, can be determined. With p2 and U2 defined 
as functions of (p^, n^, s) via figure |2] (2.12) is an algebraic equation that 
determines s as a function of (/O"*", n"*"). In general, (2.12) must be solved nu- 



merically, however, in §|5]we provide an example where an analytic solution 
is possible. Finally, we point out that (i) the restriction um > must be 
imposed, where um = s + m/ pm is the speed corresponding to the maximum 
traffic density, pM', (H) the cases of U2 = ui, where ui is defined in figure [2] 
or of no intersections between the solid and dashed curves of figure [2] do not 
yield jamitons. 

The methodology summarized above is reminiscent of the related anal- 
yses in gas dynamics [17,18], shallow water theory [19,23,24], astrophysical 
accretion flow [25], and Newtonian flow in elastic tubes [26], where ordinary 



differential equations similar to (2.3) are obtained. Indeed, (2.12) is the 



exact analog of the Chapman- Jouguet condition in detonation theory [17]. 
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Summarizing the above discussion, the following algorithm may be ap- 
plied to determine the jamiton structure: 

(i) For a prescribed downstream state, (/O^, u"*"), the regularization con- 



dition (2.121 specifies the permissible value(s) for the wave speed, s. 



(a) Once s is determined, the state upstream of the shock, [p , u ), is 



computed using the Rankine-Hugoniot conditions (2.4 - 2.5 



(in) Equation (2.3) is then integrated forward in r/, from the initial condi- 



tion n = u"*" up until u = u is reached — the traffic density, /o, follows 



automatically from (2.2). This defines the period A of the traveling 
wave or, equivalently, the circumference of the periodic roadway. 

(iv) The total number of vehicles, N, which remains fixed in time when 
there are no on-ramps or off-ramps, is then evaluated from 



[ pdx. (2.13) 

^0 



The jamitons thus obtained have the following properties: (i) The traffic 
speed smoothly increases in the downstream direction (du/drj > 0), ex- 
cept at the location of the shock, across which there is an abrupt drop in 
u; (a) The traffic density smoothly decreases in the downstream direction 
[dp/dr] < 0), except at the location of the shock, across which there is an 
abrupt increase in p. This is consistent with the experimental observations 
of Sugiyama et al. [15]. 

The above algorithm provides a parameterization of the periodic jami- 
tons using (/>"*", «"*"). Equivalent parameterizations, in terms of {p~ , u~), 
are just as easy to produce. However, these are not necessarily ideal param- 
eterizations. For example, in order to predict what jamiton configuration 
might arise from a given set of initial conditions, on a given closed roadwayj^ 
a parameterization in terms of the roadway length. A, and the total number 
of vehicles, M, would be more desirable. On the other hand, for the example 
considered in § |5] and all the other case studies that we have examined to 
date, p~ maps in a one-to-one fashion to J\f. Thus, by applying the above 
algorithm, one can iteratively determine the unique traveling wave solution 
corresponding to particular choices for the functions p and u, the parameters 
r and A, and the average traffic density Pavg 



Say to predict the patterns arising in the experiments by Sugiyama et al. [15]. 
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3 Alternative description of the momentum equa- 
tion 



The existence of self-sustained shock solutions is not specific to the details of 
the continuum model used. It is, in fact, a feature of models involving hyper- 
bolic conservation laws with forcing terms under rather generic conditions. 
To illustrate this point, we shall briefly consider the equations presented by 
Aw & Rascle [11] and Helbing [12]. These models were developed in response 
to the criticisms of Daganzo [22], who argued that second order models are 
necessarily flawed — they predict, for example, shocks overtaking individual 
vehicles in unsteady traffic flow, and negative vehicle speeds at the end of a 
stopped queue. 

Aw & Rascle [11] overcome such limitations by applying a convective, 
rather than a spatial, derivative when modeling the effects of preventive 
driving (the anticipation term in their nomenclature). This leads to the 



momentum equation (1.2) being replaced by 



{u + p)t + u {u + p)x = - (u-u) . 



(3.1) 



By substitution of the mass continuity equation (1.1 ), (3.1 ) can be re- written 
as 



ut + {u- p) Ux 



1 



{u — u) 



(3.2) 



Introducing the self-similar variable rj defined by (2.1 ), it can be shown that 

u — s) {u 



du 



u 



m 



(3.3) 



which is identical to (2.3), except that the sonic point is n ow predicted to 
occur when u — s = mP^. As before, the singularity of (3.3 1 — i.e. the sonic 
point, where the denominator vanishes — is regularized by aligning the sonic 
point with a root oiu — u. The remainder of the analysis is entirely similar 



to that outlined previously for (1.1 - 1.2), with appropriate modifications 



to the shock conditions. Under the Aw & Rascle formulation, the Rankine- 
Hugoniot condition corresponding to the conservation of momentum takes 
the form 

s\p{u^-p)\ = \pu-p^ p'u^\, (3.4) 



which replaces (2.5). 



Helbing [12] generalized Aw &: Rascle's model one step further by defin- 
ing two traffic pressures, both functions of p and u, such that the momentum 
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equation reads 

ut + uux^ ^PxH :r-Un, = -{u-u). (3.5) 

P op p OU T 

Proceeding as above, the following familiar expression can be readily ob- 
tained: 

Au {u — s)(u — u) 



where 
in which 



dry [u — sj^ — 

<^'^ = cj + mcl, (3.7) 



dpi 2 1 dp2 

— and C2 = — 2^- 

op p^ OU 



The Payne- Whitham and Aw &; Rascle results may be recovered from (3.6) 
by setting, respectively, pi = p, P2 = 0; and pi = 0, P2 = —up^ dp /dp. 



4 Phase plane analysis 

Generalizing the analyses of § [2] and § |3] dn/dry may be expressed as 

(n — s) {u — u) 



du 
drj 



(4.1) 



where $1 = 1 and m, respectively, for the Payne- Whitham and Aw & Rascle 
models. Further information regarding model behavior near the sonic point 
may be gleaned by introducing the phase plane variable ^, and rewriting 
(4.1) as the following pair of ordinary differential equations 

^ = {u- s) {u- u) = Fi{u), 



(4.2) 
(4.3) 



Note that the sonic point is a critical point of (4.2) and (4.3). The Jacobian, 
J, of the above pair of equations is then given by 

J 



Fi,u 

F2,u 



(4.4) 



12 



where 



Fi,u = -(n-s)|^^— ^ + l|+ii-n, (4.5) 
F2,u = 2{u-s) + - y2Ppp, (4-6) 



in which 



dn , d^p dc^ 

^''''"dp^= dp- 

Therefore, at the sonic point the eigenvalues (ri and r2) of J' are given by 



ri = and = ^ - O^/^ c = (4.7) 

When r2 > (respectively ^2 < 0), du/drj > (respectively du/dry < 0) 
at the sonic point. Because of the Lax entropy conditions described earlier, 
u should be a monotonically increasing function of rj away from a shock, 
which in turn requires 

r2 > <S=^> m\up\>il,c^. (4.8) 

Unfortunately, since ri = 0, the critical point is linearly degenerate. Thus 
a complete analysis of the solution behavior near this critical point requires 
a careful, but ultimately tangential, examination of the leading order con- 
tributions by nonlinearities. Nevertheless, an interesting observation can be 
made: 

As we illustrate by way of example m § r2 = coincides 

with the boundary wherein a constant uniform base state be- , . 

comes unstable to infinitesimal disturbances — see figure |4] 

and Appendix [Cj 

The theoretical underpinnings of this coincidence are not entirely clear. 
However, our numerical experiments show that there is a strong connec- 
tion between jamitons and instabilities: 

When a uniform traffic state is linearly unstable, the insta- , > 
bility consistently saturates into a jamiton-dominated state. 

For the more general analysis of Helbing [12] discussed at the end of § [Sj 
F2^u is given by 

F2,u = 2{U-S) + - r2?p, 

TTl ^(Uj — 
= 2{u-s) + - ^Pl,ppH P2,u-P2,up, (4.11) 

[u — sj^ m 
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where is defined by (3.7). Given the form of the Jacobian matrix J ^ 



modifying i<2,ti does not alter the eigenvalues r\ and r2- 

5 An example 
5.1 Preliminaries 

To make the ideas of the previous sections more concrete, we consider herein 
particular forms for p and u, and carefully examine the resulting range of 
solutions. As alluded to above, various expressions for "p and u have been 
proposed in the traffic literature. Consistent with the spirit of previous 
studies (e.g. [9,27]), our motivation is to select relatively simple functions so 
that the concepts of § |2] are succinctly illustrated with minimum algebraic 
investment. Thus, a Lighthill-Whitham- Richards forcing term of the form 

n = nofl-— (5.1) 
V Pm) 

is chosen. Moreover, by analogy with the shallow water equations [23], we 
select 

P=i/3p', (5.2) 

so that -pp = (? <yi p [11]. Alternatively, one could define p to be singular 
in /9, such that pp oc p'^'^ /{pM — p)'^'^ for some ^i, (p2 > 0. Whereas this, 
more complicated, choice for p enforces pM as the maximum traffic density, 
the resulting algebraic relations become somewhat unwieldy. We therefore 
defer consideration of singular pressure functions to future studies. 

Applying the above definitions to the Payne- Whitham model of §[T]yields 



dv {u - s)2 - 



(5.3) 



(u-s) 

The sonic point is then defined by 

n-s = (^m)^/^ (5.4) 

while the zeros of the numerator are given by 

1 / \ 1 / \ { 4 ttQ m 1 "^^^ 
u = {ui, U2} = s + ^{uo - s) ± i,{uo - s) <1 — ^ . (5.5) 
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The regularization condition (2.12) may then be written as 

U2 = s + {l5mY'^. 



(5.6) 



Equations (5.5) and (5.6) yield a cubic polynomial in s, with at most two 
physically-relevant roots. With p defined by (5.2), the Rankine-Hugoniot 
conditions specified in (2.4) and (2.5) may be combined to yield 

i{-l + (l + 8M^ 



PI 



Vf2)l/2| 



where M_ 



(5.7) 

— s\/{j3 p^Y^"^ is the upstream Mach number [13,14]. 



Finally, the non-trivial eigenvalue, r2, of the Jacobian matrix given by (4.4), 
is, at the sonic point. 



r2 



PM V /3 



2\ 1/3 



{I3m 



,1/3 



(5.8) 



Having regularized the ordinary differential equation given by (5.3) — 
i.e. upon defining U2 by (5.6) and canceling u — U2 from the numerator and 
denominator, the resultant differential equation has an exact, albeit implicit, 
solution given by 



V 



(/?m)^/^ , / u — s 
u + — — In ' 



Ul — s 
+ { {(3m)'^/^ + ui -s + 



— s ^ 
(^m)2/3 
Ul — s 



In 



Ul — u^ 
Ul — u 



(5.9) 



A similar result is obtained in the study of roll waves in shallow water 
systems - see e.g. (4.18) in [23]. 

By definition, u = u~ when ij = X, where A is the length of the periodic 
roadway. Therefore 



A 



u 



u~ + 



Wmfl^^ (u 
in ' 



Ul 



s 
— s 

Ul — s 



In 



Ul — u' 
Ul — u~ 



(5.10) 



Starting from (2.13), the total number of vehicles along the periodic circuit 
can be computed from 



pdx 



pdrj = m 



+ u 



di]) 



(5.11) 
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Application of (5.3) in ( 5.11[ ) yields then the explicit result 



m T 



(/?m)2/3 r u 



u ' 



Ul — s \{u — s) (li+ — s) 



+ 



til — s 1 Ul — s J — s) {ui — w) 

Ul — u'^^ 



+ In 



Ul — u 



(5.12) 



Comparable exact solutions may also be determined when p oc p , rather 
than p (X p'^. These are presented in Appendix [Xj 



5.2 Numerical method 

In order to validate the aforementioned theoretical solutions and assess j ami- 
ton stability, we performed numerical simulations using a Lagrangian par- 
ticle method [28]. In this method, each discrete particle, i, is assigned an 
initial position, Xj, and speed, Ui, which subsequently changes in time ac- 



cording to (5.13 - 5.18). The mass balance equation (1.1 ) is satisfied identi- 



cally as the particles move, i.e. the numerical scheme is mass-conservative by 
construction. Importantly, the number of particles is typically two orders of 
magnitude larger than the number of vehicles. Thus, although a Lagrangian 
approach is employed, the numerical scheme constitutes a macroscopic, not 
a microscopic, description of traffic flow, albeit one with an intuitive link 
between the particle and vehicle density. 

In general terms, the numerical method solves the differential equations 



(5.13) 




where the particle acceleration, -iij, is expressed as 
du . 



dt 



{Xi, t) = Ut{Xi, t) +UUx{Xi, t), 



and 



Oi = 7 7^ Px{Xi, t) + - {u{p{xi, t)) - u{Xi, t)} . 

p{Xi, t) T 



(5.14) 



(5.15) 
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The distance between adjacent particles i and i + 1 is defined by 
rEj+i — Xi. Then the inter-particle density is computed by 



2 



where = N /Np — in which J\f is the number of vehicles (as specified by 



(2.13)) and Np is the number of particles. From (5.16), we define the vehicle 
density and the density gradient using a non-equidistant finite-difference 
stencil [29]: 



1 



(5.17) 



^^'"^ min{d._^i, 2a!^_i} +min{(i._i, 2d._^i}' ^^'^^^ 

The denominator in (5.18) is chosen so that, at the location of the shock, 
a given particle is influenced only by its nearest neighbor. The numerical 
scheme is stable — however, it produces bounded but unphysical oscillations 
wherever the vehicle density or speed changes abruptly. We suppress these 
features by adding a small amount of numerical viscosity, so that at each 
time step the velocity profile is smoothed out. 

5.3 Results 

Here we compare the theory of self-sustained traffic jams developed in §[2]and 



5.1 and numerical solutions of ( |l.l| - [r!2| ) — obtained using the algorithm 



of § 5.2 Employing the forcing and traffic pressure terms specified in (5.1) 



and (5.2), the equations have the non-dimensional form 

p*t* + {u* p*)^* = 0, 
Ti «+tx*<,) + (r2/p*)p*. = l-p*-n 



(5.19) 



where Ti = tuqpm, ^2 = P'Tp\i/uQ^ and p* = \{p*)'^. The non-dimen- 
sional and dimensional variables are related via 

P = Pm P* 1 u = uou*, x = £x*, and t = {£/uo)t*, (5.20) 

where £ = 1/pm is the effective vehicle length, defined earlier in § [T] the 
non-dimensional traveling wave speed is s* = s/uq, and the parameter r 
is the adjustment time-scale required for an individual vehicle to effect an 
0(1) change in speed. 
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We consider as canonical modeling parameters i = 5m, uq = 30m/s, 
r = 3.3s, and (3 = 450m^/s'^, whereby 

Ti = 20 and T2 = 2. (5.21) 

Note that (3 = / p is related to the speed at which disturbances propagate 
through traffic. Because a description, both theoretical and numerical, of 
jamitons is the central focus of the present analysis, we deliberately select a 
value for /? towards the lower end of its representative range, so as to facili- 
tate an explicit display of jamiton properties and behavior. Hence, /?, whose 
exact numerical value is likely difficult to estimate in any event, is chosen so 
that infinitesimal perturbations to a uniform base state develop into jami- 
tons when the base-state density is 10% of pM — see Appendix IC] and. 



more particularly, the stability condition (5.22). Whereas larger numerical 
values for /3 could be selected, corresponding to a more restrictive insta- 
bility condition, the jamitons observed in these cases, though qualitatively 
equivalent to those described below, are somewhat less easy to visualize. 
The practical ramifications associated with choosing a larger value for /3 are 
briefly addressed when we discuss (after examining the results in figure [4| 
the possibility of vehicular collisions. 

The steady-state variations of u* and p* , as functions of the non-dimen- 
sional variable r] r pM = x* — s* t* , are shown in figure |3] for a circular road 
of length A = 100^, with two different choices for the conserved number of 
vehicles, M. The shock occurs at the two extreme ends of the horizontal 
domain, and connects the ratios /uq to u'^/uq, and p^ / pM to p'^/pM- In 
figure |3]a, M = 27, and the maximum traffic density (i.e. p~^) is predicted 
to be just below pM- Conversely, in figure |3]b, M is increased to 38, which 
results in > pM- The physical implications of this result are considered 
in the following two paragraphs. Both theoretical and numerical data are 
included in figure [3j The comparison is very favorable, except right at the 
shock location — where the numerical scheme smears the shock. 

Figure |4] indicates, as a function of the normalized average traffic density 
Pavg/PM (where pavg = AT/A), the range of possible solutions allowed by the 
model equations in § 5.1 The thick solid curves of figures [4]a,b show p~ j pu 



and p'^ / pmi respectively, whereas the thick dashed curves show u~ /uq and 
u'^/uq, respectively. As with the discussion in § [2| we observe that, for 
a prescribed road length A and model parameters pM, uq, /?, and r, the 
average density pavg (hence M) uniquely determines the flow conditions to 
either side of the shock. Although the solid and dashed curves of figure |4]a 
do not proffer any particularly meaningful insights, those of figure |4]b are 
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Figure 3: Comparison of theoretical (thick curves) and numerical (thin 
curves) solutions. The numerical solutions are the final, asymptotic state 
of an evolution started with a small perturbation of a uniform unsta- 
ble base state — see (5.22). The agreement of this final state with the 
jamiton corresponding to the same road length, A, and number of vehi- 
cles, M, is remarkable. The plots show profiles of p/pM and u/uq versus 
V' Pm = {x — st) pm- The equations and parameters are as in ( 5.19| - 5.21), 
with A = 100^ = 500 m. In panel (a) N = 27 {pavg/pM = 0.272). In panel 
(b) AA = 38 (/5a^g/pA/ = 0.384). 



significant in that they predict the following forms of model breakdown: 
P^ > Pm when pavg/pM > 0.277 and < when pavg/pM > 0.391 

Rather than signifying a fundamental modeling flaw, observations of 
model breakdown provide helpful guidance in the safe design of modern 
roadways. In particular, for the choice of parameters germane to figure |4j 
vehicular collisions (i.e. p — > pm) are anticipated once the average center- 
to-center separation between adjacent vehicles falls below approximately 
3.6^. Needless to say, this simple result is not universal for all types of 
traffic flow or roadway conditions. We consider herein a periodic track of 
prescribed length (500 m), particular functional forms for the traffic pressure, 
p, and equilibrium speed, u, and a liberal numerical value for /3 such that 
jamitons appear even in relatively light traffic. In particular, choosing a 
larger value for /? would delay, though not necessarily avoid, the onset of 
vehicular collisions. 

The roadway's carrying capacity could, in principle, be increased if two 
or more traveling waves were to be forced rather than the single density spike 
considered here. However, our numerical simulations suggest that such bi- 

^In Appendix B we verify that p — > pM before u"*" 0. 
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Figure 4: Parameters characterizing the exact jamiton solutions in 
a fixed road length A = 500 m, as functions of the number of vehicles, N . 
Specifically, curves of p~ I pM-, I Phi-, u" /uq, and u^/no versus Pavgl Pm = 
M l{X Pm) are shown, in addition to r2/tto — as given in (5.8). The vertical 



dashed line is the stability boundary specified by (5.22). The equations and 



parameters are as in (5.19-5.21). Note that (i) r2 goes through precisely 



at the stability boundary — see (4.9); (ii) the peak density p"*" reaches pu 
before u"*" becomes negative — see Appendix pi 



or multi-modal structures are often unstable, and that they quickly coalesce 
into a single structure. By contrast, unimodal traveling wave solutions of 
the type illustrated in figure [3] appear to be stable to small perturbations — 
except in the special case where the periodic roadway is made to be quite 
long, i.e. ~ 2 km for the parameters appropriate to figures [s] and |4 On an 
extended circuit, the traveling wave solutions have a long, nearly constant 
state downstream of the shock — whose average density exceeds the require- 
ment for linear stability discussed below. Thus infinitesimal perturbations 
may grow, leading to further traveling waves, albeit of slightly different 
amplitude from the original. Owing to the circuit length (i.e. ~ 2 km, as 
compared to the 500 m long circuit considered in figure [s]), these traveling 
waves are dynamically independent in that wave coalescence does not occur 
for a rather long time, if at all. 

While the above observations would benefit from the development of 
a more fundamental framework, it can be difficult to describe analytically 
the spatio-temporal stability of self-sustained traveling waves, due primarily 
to the transonic nature of the solution (see e.g. Stewart & Kasimov [18], 
Balmforth & Mandre [24] and Yu & Kevorkian [30]). Work is on-going 
to adapt nonlinear stability analyses for specific application to continuum 
traffic models. 
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In Appendix [C] we compute the boundary between linearly stable and 
unstable uniform base states, which is given by 



P 

PM 



Ppm 



(5.22) 



In connection with this boundary, we point out the following facts, which to- 



gether with (4.9-4.10), reinforce the point that there is a strong connection 



between jamitons and uniform flow instabilities. 

• Figure [4] shows that [p] ^ and [u] ^0, as Pavgl Pm i Pm/u^- Thus 
the jamiton amplitude becomes smaller as the corresponding uniform 
flow becomes less unstable, and vanishes altogether in the limit. 



The stability criterion (5.22) also specifies the location at which ui 



U2, where ui and U2 are defined in figure [2] Once ui and U2 coalesce, 
no self-sustained shock wave may occur — since condition (a) in (2.11 ) 



fails. This is consistent with the discussion of the previous bullet, of 
course. 

Thus jamitons become possible when the corresponding uniform state be- 
comes unstable and — see (4.10 ) — the basic flow state changes from uniform 



flow to a jamiton-dominated state. In other words, evidence indicates that 



A crucial bifurcation in the traffic flow behavior occurs at the 



stability boundary prescribed by (5.22). 



(5.23) 



We defer a more in-depth investigation of this question for future work. 



Finally, note that although the boundary specified by (5.22) is mathe- 



matically robust, in practice it may become "fuzzy" owing to the possible 
breakdowns of the continuum hypothesis, especially at low vehicle concen- 
trations. 



6 Conclusions 

In this work, we have found and successfully exploited a strong similarity 
between gas-dynamical detonation waves and shocks in traffic flow, in order 
to develop a theory of "jamitons," steady self-sustained traffic shocks. Jami- 
tons naturally arise from small instabilities in relatively dense traffic flow, 
and can be interpreted as saturated phantom jams. While a single jamiton 
may not necessarily significantly delay individual vehicles, a succession of 
jamitons, as might arise during rush hour, for example, is expected to frus- 
trate motorists over long lengths of roadway. Moreover, jamitons represent 
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regions in which the traffic density increases dramatically over a relatively 
short distance [15] and, as such, are hot spots for vehicular collisions. 

The analogy drawn above goes beyond the weaker analogy with inert 
shock waves in gas dynamics considered by Kerner, Klenov &: Konhauser 
[27]. Unlike inert shocks, detonation waves can be self-sustained, due to the 
existence of a sonic point in direct correspondence with the traffic model we 
consider here. 

Using the Lax entropy conditions for hyperbolic conservation laws, we 
show that for realistic pressure and equilibrium-speed functions, the only 
allowable self-sustained shocks are those that are overtaken by individual 
vehicles. Moreover, for the simple, but widely-applied, choices u = uo{l — 
p/pm) and p oa p'^, with 7 = 1, 2, we are able to describe the jamiton 
structure analytically. Theoretical solutions show excellent agreement with 
the output from direct numerical simulations of the governing system. 

Examples of jamitons with 7 = 2 are illustrated in figure |3j Beyond a 
critical average density, we predict solutions for which the density immedi- 
ately downstream of the shock exceeds the maximum allowable density, pM, 
corresponding, physically, to a state of vehicular collisions (figure |4|. Such 
instances offer important design insights for roadways; most obviously, it is 
advantageous to choose speed limits and roadway carrying capacities so as 
to avoid circumstances where densities with p ~ pM are "triggered" (say, by 
a jamiton) anywhere within the domain. 

Having identified self-sustained traveling wave solutions in traffic flow, a 
major objective of future research is to ascertain their spatio-temporal sta- 
bility. Some of the analytical challenges associated with this line of inquiry 
are identified in § |5.3[ Resolving issues of stability offers the possibility 
of increasing roadway efficiency, for example, by exciting multiple traveling 
waves of relatively low density as compared to the single density spike ex- 
hibited in figure |3j We hope to report on the results of such an analysis 
soon. 
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A Exact solutions when p (x p 



Equations (5.9), (5.10) and (5.12) specify, respectively, exact analytical so- 
lutions for n, A and J\f when p (x p^. Corresponding expressions, valid when 
p (X p, are given by 

rj = u'^ - u+ (ui- s + In 
X = u+ -u- + (^ui- s + In 



Ul — 
Ul — u 

Ul — li^ 

Ul — w 



and 
M = mr 



+ 1 + 



In 



Ul — u' 



(A.l) 
(A.2) 

(A.3) 



Ml — S \U^ — S J \ Ul — S / \Ui — U 

respectively. As before, the traffic speed is given as an implicit function of 



r] in (A.l ). 



B Model behavior for p > pM 

From the discussion of § [sj vehicular collisions are forecast once = pM, 
whereas nonsensical vehicle speeds are predicted once m"*" < 0. It is demon- 
strated herein that the former condition is necessarily achieved before the 
latter. 



Equation (2.2) shows that 



PM {UM 



m > 0, 



(B.l) 



where um, which must be positive for jamitons to exist, is defined in figure 
[2]b. As M^ 0, the left-hand side of (B.l ) approaches —p^ s, demonstrating 
that s < 0. In this limit, therefore. 



EL 

PM 



> 1. 



(B.2) 



This result can be understood intuitively by examining the functional form 



of the Lighthill-Whitham- Richards forcing term (see (5.1)): n < requires 



p > pm- Clearly, the sensible alternative is to define u = ior p > pM- The 
point is moot, however: this amounts to correcting the model equations in 
a regime that is already physically unrealistic. 
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C Linear stability of the Payne- Whit ham 
model considered in § [5] 



For a right-hand side forcing function of type (5.1), the constant base state 

(C.l) 



solution to (1.1) and (1.2) is given by 

u = uq I 1 



P = P, 



P 

PM 



Following ideas discussed in Kerner &; Konhauser [9], Helbing [12] and else- 
where, the linear stability of this base state can be explored by introducing 
perturbation (hatted) quantities, defined such that 



P 



U = Uq (1 '■ — ) -|- U, 

PM, 



(C.2) 



P = P + P, 

where p and u are expressed in terms of normal modes by 

Here k is the horizontal wave number and a is the corresponding growth 
rate. Application of (C.2) and (C.3) into (1.1) and (1.2) shows that 



^0 , . a 



rpM 



i k p 




' R ' 




' " 


a + ik'ip + ^ 




U 








(C.4) 



where, for notational economy, we have introduced 

ip = uq(i — 

V PM 



(C.5) 



Requiring that the determinant of the matrix from (C.4) vanish shows that 

1 



a 



-ikij- — il + T), 

2 T 



in which 



= l-Ak^PpT'^+A\kuQT 



PM 

Generically, T may be written as T = Ai — iA2, where 

A? - Ai = 1 -4fe2/3pr^ 



(C.6) 
(C.7) 

(C.8) 
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and 



A1A2 = -2A:uoT— . (C.9) 

PM 



Eliminating A2 from (C.8) and (C.9 1 yields the following polynomial expres- 
sion: 

P(A2) =Af-(l-4A:2/3pr2)A2-4A;2^t2^2^ = 0. (C.IO) 

Pm 

Linear stability requires a non-positive growth rate, i.e. 

Real((T) < < A^ < 1 V{1) > 0. (C.ll) 



From (C.IO), the latter condition is satisfied provided 



/3>^. (C.12) 
Pm 



This completes the derivation of (5.22). 
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